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I. INTRODUCTION 


Single and few-layers of transition metal dichalcogenides show promising electronic, op¬ 
tical, and mechanical properties and are considered as potential candidates for future elec¬ 
tronic applications. 20 Because of weak inter-layer van der Waals bonds in their layered struc¬ 
ture, single to few-layers of these materials can be obtained by mechanical or chemical 
exfoliation techniques. 2,21,27 Among these materials single-layer M 0 S 2 has attracted the at¬ 
tention of scientists. 13,18,19,29 Single-layer M 0 S 2 has a direct band gap of 1.8 to 1.9 eV, 19,29 
which makes it suitable for various electronic applications. 32 It has been shown that the 
application of compressive and tensile biaxial strain results in an indirect bandgap. 0,8,30 An 
MOS transistor based on this material has demonstrated a I on /I Q s ratio of 10 8 , a relatively 
steep sub-threshold swing of 74 mV/dec and an extremely small off-current of 25 fA//im, 20 
moreover possible applications to hetero-junction inter-layer tunneling FETs have also been 
proposed and theoretically investigated. 15 Room temperature mobility of n-type single-layer 
MoS 2 has been reported to be in the range of 0.5-3 cm 2 /(Vs) and can be increased to about 
200 cm 2 /(Vs) with the use of high-ft dielectrics. 26 

The low-held mobility is one of the most important transport properties for a large 
number of physical systems and electronic devices. A comprehensive study of strain effects 
on the mobility of single layer M 0 S 2 , however, is missing. I 11 the present work, the effects of 
biaxial and uniaxial strain on the low-held mobility of single-layer M 0 S 2 is investigated by 
using ab-initio simulations along with the linearized Boltzmann transport equation (BTE). 23 
Scattering rates due to intrinsic phonon, charge impurities, and remote phonon are taken 
into account. 


II. BANDSTRUCTURE AND SCATTERING RATES 

In the hrst part of this section, some details about the ab-initio calculations of the elec¬ 
tronic bandstructure in the presence of strain are discussed. Thereafter, the formulation of 
various scattering rates is described. 
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A. Bandstructure 


We carry out first-principle simulations based on the density-functional theory (DFT) 
along with the local density approximation (LDA) as implemented in the SIESTA code 9,14,28 
to investigate the relevant electronic properties of a single layer M 0 S 2 under strain. While 
DFT-LDA, in general, underestimates band gaps, the resulting dispersion of individual 
bands, i.e., effective masses and energy differences between valleys, is less problematic. 10 A 
cutoff energy equal to 600 Ry is used and a vacuum separation of 30 A is adopted, which is 
sufficient to hinder interactions between adjacent layers. Sampling of the reciprocal space 
Brillouin zone (BZ) is performed by a Monkhorst-Pack grid of 18 x 18 x 1 h-points. Calcu¬ 
lations begin with the determination of the optimized geometry, that is the configuration in 
which the residual Hcllmann-Feynman forces acting on atoms are smaller than 0.01 eV/A. 
The calculated lattice constant of unstrained single-layer M 0 S 2 is 3.11 A that has a good 
agreement with the reported value in Ref. 1 and 3. Fig. 1(a) shows the energy contours 
of the conduction-band in the Erst BZ for unstrained single layer M 0 S 2 . In an unstrained 
material the lowest and the second lowest minimum in the conduction band are denoted as 
K-valley and Q-valley, respectively. The 6 K valleys are degenerate in the unstrained and 
in all the strained conditions explored in this paper. The 6 Q-valleys are degenerate in un¬ 
strained conditions, while under uniaxial strain they split into 4 Qa- valleys and 2 QB-valleys 
with different effective masses and energy minima, as discussed and illustrated in Sec. IV. 
The energy distance between K-valley and Q-valley for unstrained material is evaluated to 
be 160 meV, in agreement with Ref. 10. Fig. 1(b) shows the calculated DFT-LDA band 
structure and depicts a direct band gap of 1.92 eV at the K-point which is very close to the 
experimentally measured value of 1.85 eV. 29 


B. Scattering with M0S2 phonon modes 

Scattering rates due to intrinsic phonons (including acoustic, optical and polar-optical 
phonons), to remote phonons and to charged impurities are taken into account. Piezoelectric 
coupling to the acoustic phonons is only important at low temperatures and is neglected 
in this work. 11 If the surrounding dielectric provides a large energy barrier for confining 
electrons in the M 0 S 2 layer, the envelope function of mobile electrons can be approximated 
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as 4/ k{r,z ) = x( z ) exp (ik.f^j/y/S with x( z ) — \j (2/a) sin(7rz/a), 17 where S' is the area 
normalization factor, fc is the in-plane two-dimensional wave vector, a is the thickness of 
single layer M 0 S 2 and r is the in-plane position vector. The scattering rates for the acoustic 
and optical phonon are discussed first. 

Using Fermi’s golden rule the scattering rate from an initial state k in valley v to the 
final state k' in valley w can be written as 

_ _ 2tt - -> 

S v ’ w (k, k') = —\M v ’ w (k, k')\ 2 8[E u \k') - E v (k ) T hu(q)\ , (1) 

where \M v,w (k, k')\ is the matrix element for the mentioned transition and hoj(q) is the 
phonon energy that may depend on q = \k — k'\. The intra-valley transitions (v = w) 
assisted by acoustic phonons can be approximated as clastic and the rate is given by 

S,c(£, k') = - E(k )] , (2) 

where k b is the Boltzmann constant, T is the absolute temperature, _D ac is the acoustic 
the deformation potential, p = 3.1 x 1CU' [gr/cm 2 ] is the mass density and v s is the sound 
velocity of single layer M 0 S 2 . On the other hand, the rate of inelastic phonon scattering, 
including intra and inter-valley optical phonons, and inter-valley acoustic phonons, can be 
expressed as 


tt(D v,w l 2 

qv,w (t Cr\ _ V ac/op/ 
°ac/opl K ’ K ~ ,, n Q 

^ac/op 


^-op ^ a "F n 


5{E w (k') - E v (k) T h(u ac/op (g)] , 


(3) 


where D^J op is the acoustic/optical deformation potential for a transition between valleys v 
and w, hco ac / op (q ) is the phonon energy, and n op is the phonon occupation (upper and lower 
sign denote phonon absorption and phonon emission, respectively). The phonon assisted 
inter-valley transitions considered in this work, and the corresponding phonon momentum 
are shown in Fig. 2. In our calculations, we employed the deformation potentials and phonon 
energies from Ref. 16 that are reported in Table I and Table II. It should be noted that the 
same deformation potentials are used for Qa and Qb valleys. 


C. Remote Phonon Scattering 

Another important scattering source considered in this work is the remote phonon or 
surface-optical (SO) phonon scattering mechanism. The source of this scattering is in the 
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surrounding dielectrics via long-range Coulomb interactions, provided that the dielectrics 
support polar vibrational modes. By assuming semi-infinite oxides and neglecting the pos¬ 
sible coupling to the plasmons of the two-dimensional material, the energy dispersion of SO 
phonons can be obtained by solving the secular equation 22 


( e box(^) + OdX^oxXX + Od) — XboxtX) — OdX^oxXX — OdX ~ 0 j 


(4) 


where the thickness a of the single layer MoS 2 is set to a = 3.17 A , 33 c 2 d is the dielectric 
constant of the two-dimensional material (single layer MoS 2 in this work), the index box 
and tox denote the back-oxide (z < 0) and the top-oxide (z > a), respectively, q — \k — k'\ 
and e 2 D for MoS 2 is set to 7.6. 17 A numerical solution of Eq. 4 shows that the frequency 
of remote phonon has a very weak dependence on q, that consequently we neglected in our 
calculations by setting e~ 2qa ~ 1 in Eq. 4. With this approximation, Eq. 4 simplifies to 
ebox(^) + etox(u>) = 0, that we solved by using the single polar phonon expression for the 
e ox (a;) in each oxide: 

,0 _ ,oo 

e °*M = e°° + 2/2 ’ ( 5 ) 

-L UJ j UyrpQ 

where e°° and e° are the high and low frequency dielectric constant, respectively, and (Uto 
is the frequency of the polar phonon in the oxide. We could provide analytical solution 
for Eq. 5 and express cu so , b ox as: uj 2 q box = (-B + \JB 2 - 4AC)/(2A) and for ^so,tox as 


Wo,tox = (-B - VB 2 -4AC)/(2A), where A = (e~ x + eggj, B = -(e° ox + eggjw, 


'so 


2 

TO,tox 


( e box + e tox)^TO,box and C = (C + ebox) w TO,tox w TO,box- Table 111 reports the parameters of 
dielectric materials that are studied in this work and indicates the corresponding calculated 
SO phonon frequencies. The scattering matrix element of remote phonon can be written 


as 


.22 


M SOjtox (fc, k f ) = 

■^SO,box(^5 ^ ) 


^^SO,tOX j 

( 1 

1 \ 

2 Sq 1 

^ ^tox ^box(^so,tox) 

^tox ^box(^so,tox) / 

^^sOjbox j 

( 1 

1 ) 

2 Sq 1 

^ ^box ^tox(^so,box) 

^box ^tox(^so,box) / 


( 6 ) 

(7) 


Scattering with SO phonon mode is inelastic and we consider only intra-valley transitions. 
The corresponding transition rate is 


2vr, 


S so {k ) k') = ^\M so {k,k')\ 2 


1 1 
n so + 2^2 


8[E(k') — E(k) =F huj so ] , (8) 

where n so and Ticu so are the SO phonon occupation number and energy, respectively. 
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D. Scattering with Coulomb Centers 


To investigate the effect of the dielectric environment on the scattering of carriers from 
charged impurities located inside the single layer M 0 S 2 , we assume that the charged im¬ 
purities are located in the center of the single layer M 0 S 2 thickness, that is at z = a/2. 
The Fourier transform of the scattering potential due to a charged impurity located at 
(r, z) = (0, a/2) can be written as 5 

p 2 

<j>{q, z ) = - [ e -«l*-«/21 + Ce qz + De~ qz ] , (9) 

2^620 


where e is the elementary charge, and C and D can be written as 

( e 2D - ebox)( e 2D - etox ) e_<?a/2 + ( e 2D + e° ox )(e 2 D - £to X ) eaq/2 - ( e box _ e 2D)(et ox + 62d) 


c = 


and 


( e 2D + e box) ( e 2D + £tox) e2aq - ( e 2D “ «w)(e 2 D ~ <Hox) 

(^ ox - e 2D ) [c + e-W2] 


D = 


“box 


+ ^2D 


( 10 ) 


(ii) 


Using Eq. 9, the x( z ) form x( z ) — \J (2/a) sin(7rz/a) and assuming intra-valley transitions 
for scattering with charged impurities, the transition matrix elements take the form 
.2 


M^(k,k') 


q 


qae 2 D \q q 2 + (27t /a) z 

+ ^(l + " 


x 


x (<"’ -1) + k (1 _»'■) _ 


( 12 ) 


qae 2 o \q q 2 + (2 k/ a)' 

where q = \k — k'\. Eq. 12 expresses the matrix element for a Coulomb center located in 
(r, z) = (0, a/2) and does not account for the screening produced by the free carriers in MoS 2 ; 
such a screening effect is introduced according to the dielectric matrix approach discussed 
in Sec. HE. The overall matrix element produced by a set of Coulomb centers randomly 
distributed at positions (r, a/2) is known to be affected by the statistical properties of the 
distribution and, in particular, by a possible correlation between the position of Coulomb 
centers. In this paper we do not address these difficulties and simply write the overall matrix 
element as | M cb (k, k ')| J = N D \M^\k, k')\ 2 /S, where N D is the impurity density per unit 
area and is given by Eq. 12. Scattering charged impurities is treated as elastic and the 
rate is therefore given by 


j ^ / 'j / _^ ^ ^ j 

Scb(M') = -j-\M ch {k,k')\ 2 5{E{k') -E{k)) . 


(13) 


6 














E. Screening 


The effect of static screening produced by the electrons in the MoS 2 conduction band is 
described by using the dielectric function approach, 5 so that the screened matrix element 
M™ T (k, k') in valley w is obtained by solving the linear problem: 


w 


(14) 


where M v (q) are the unscreened matrix element. As can be seen in Fig. 1(a), there are three 
different valleys in strained single layer MoS 2 (K, and Q b valleys), hence v, w £ {K, 
Qa, Qb}- In Eq. 14, e v,w is the dielectric matrix which is introduced as: 

e V ’ w (q ) = 5,,. - ^ -rn w (q)F v ’ w (q) , (15) 

q(e 2D + e box ) 

where S lKW is the Kronecker symbol (1 if v = w, otherwise zero), IT"(g) and F v,w (q) are the 
polarization factor and unit-less screening form factor, respectively. 5 In the case at study 
the dielectric matrix can be analytically inverted to evaluate screened matrix elements as: 


AC(«) 


(l - E„/„«”'”’(<()) M”(q) + E„ # „ «”'”(« )M"{q) 

2-E 


(16) 


The static dielectric function approach described above has been directly used for the 
scattering due to charged impurities, while the situation is admittedly more complicated 
for phonon scattering. For the inelastic, intervalley phonon transitions described in Ta¬ 
ble I and Table II the relatively large phonon wave-vector (see also Fig. 2) and the non-null 
phonon energies suggest that it is safe to leave these transitions unscreened, because the 
dynamic descreening and the large phonon wave-vectors make the screening very ineffective. 
Arguments concerning screening for intra-valley acoustic phonons are more subtle and con¬ 
troversial and a thorough discussion for inversion layer systems can be found in Ref. 7. We 
here decided to leave also intra-valley acoustic phonons unscreened, which is the choice em¬ 
ployed in essentially all the studies concerning transport in inversion layers that the authors 
are aware of. The screening of the SO phonon scattering is also a delicate subject, because 
the polar phonon modes of the high -a dielectrics can couple with the collective excitations of 
the electrons in the MoS 2 layer and thus produce coupled phonon-plasmon modes, 7,22 whose 
treatment is further complicated by the possible occurrence of Landau damping. 7,22,31 In this 
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paper we do not attempt a full treatment of the coupled phonon-plasmon modes, 7,22 but in¬ 
stead show in in Sec. IV results for the two extreme cases of either unscreened SO phonons 
or SO phonons screened according to the static dielectric function. We can anticipate that 
while the inclusion of static screening in SO phonons implies a significant mobility enhance¬ 
ment compared to the unscreened case, the mobility dependence on the strain and on the 
dielectric constant of the high-K dielectrics is not significantly affected by the treatment of 
screening for SO phonons. 


III. MOBILITY CALCULATION 

Acoustic, optical, polar-optical, remote phonon, and charged impurity scatterings are 
considered for the calculation of low-held mobility. As it will be discussed in the next section, 
the bandstructure for and Qb valleys is not isotropic and the mobility shows direction- 
dependence, hence we calculated mobility by solving numerically the linearized Boltzmann 
Transport Equation (BTE) according to the approach described in Ref. 23, which does not 
introduce any simplifying assumption in the BTE solution. In particular, mobility has been 
calculated along the armchair and zigzag directions and strain has been also studied for the 
uniaxial configuration along either armchair or zigzag direction, as well as for the biaxial 
configuration. 

In order to describe in more detail the mobility calculation procedure, we first recall that 
the longitudinal direction of Q^-valley is neither the armchair nor the zigzag direction, and 
Fig. 1(a) shows that 9 is the angle describing the valley orientation with respect to the zigzag 
direction in k-space (i.l. armchair direction in real space). Let us now consider Erst the case 
of the mobility of valley v along the armchair direction, that can be written by definition 
as = J^/Fa, where is the current component in the armchair direction for the 
valley v induced by the electric held Fa along armchair direction. The current can be 
expressed as J^ !> = cos ( 9 V ) + J sin ( 9 V ) in terms of the current components j\ v \ J ( t v} 

along, respectively, the longitudinal and transverse direction of the valley v. By denoting the 
longitudinal (T)) and transverse component ( F t ) of the electric held as Fi = Facos[9 v ) and 
F t = FAsin(9 v ), the currents jj' v> and in turn can be written as j\ l) = Fj + F t 
and F[ + F t: where fin, /% and /% are the entries of the two by two mobility 


matrix in the valley coordinate system. 5 Consequently we finally obtain: 


(u) 


J 


(v) 


= iff cos 2 (9 V ) + sin 2 (9 V ) + sin (9 V ) cos ( 9 V ) . 


» 


(^) 


F. 


a 


(V) 


By following a similar procedure, the mobility iff of the valley v along the zigzag direction 
can be written as: 


( V ) 

/4 = 


J 


( v ) 


= ff sin 2 (9 V ) + ff’ cos 2 (0„) — 2/i^ ; sin (9 V ) cos (9 V ) . 


» 


(«) 




(18) 


For the circular and elliptical bands employed in our calculations (see Fig. 3(g)), iff is zero 
for symmetry reasons. 5 After calculating the mobility for each valley, the overall mobilities 
fj,A and fiz are obtained as the average of the mobility in the different valleys weighted by 
the the corresponding electron density. 

Eq. 17 and Eq. 18 allow us to calculate the mobility iff and iff from the longitudinal 
//]) and the transverse mobility /i)) of the valley v which are the mobilities obtained from 
the linearized BTE when the electric field is either in the longitudinal or in the transverse 
direction of the valley v. As already said, the fi v u and n\ t have been obtained by using the 
approach of Ref. 23, whose derivation for the case at study in this work can be summarized 
as follows. The out of equilibrium occupation function f v (k ) for the valley v in the presence 
of a field F x is written as 

nk) = fo(E v (k))~eF x g\k) , (19) 


where / 0 (F) is the equilibrium Fermi-Dirac distribution function, x € {Z,i} is either the 
longitudinal or the transverse direction of the valley and k = (. ki , k t ) is the wavevector in 
the valley coordinate system. Eq. 19 is a definition of g v (k), which is the unknown function 
of the linearized BTE problem. For a two-dimensional system the linearized BTE can be 
written as 23 

g v (k) ( 2 ^ J] [ A v ’ w (k, k')S[E w (k') - E v {k) T ticu ph }dk'\ 

\ n w k ' J (20) 

f e A ” - ”A k')g“(k')S{E"-(k') - E v (k) =F M dk' = vl(k) , 


where v x is the x component of the group velocity of valley v and, for convenience of notation, 
we have introduced the quantity 


A v ' w (k,k') 


\M v,w (k, k')\ 2 


1 - f 0 (E w (k')) 
1 - f 0 (E-(k)) 


( 21 ) 
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To numerically solve Eq. 20, we employed the discretization scheme introduced in Ref. 23: 
k is discretized according to a uniform angular step A/3 and also a uniform energy step. 
The discrete values k V} r td of the wave-vector magnitude correspond to one of the discrete 
energy values and the generic discrete wave-vector k Vi1 - id = (k Vjr , d , dA/3) is identified by the 
magnitude k V}Tt d and the angle dA/3 (with d being a positive integer number). For each 
scattering mechanism, by converting the integral over k in an integral over the energy and 
the angle [3 and then using the above mentioned discretization, Eq. 20 can be rewritten as: 


9 i.kv,r,o 


Al 

27 rh 


E A w,r',d' rw,r' 
^v,r,d u v,r,d 


w,r' ,d' 


Al 

2irh 


£ BgfAKsSK 


r',d' 

r.d 


w,r' ,d' 


Vx(ky,r,d)f( E ( k v,r,d))[i- ~ f(E(kv t r,d))] 
knT 


( 22 ) 


Eq. 22 is a linear problem for the discretized unknown values g(k v ^d) written in terms of 
the coefficients A^^f’ and EA'lf defined as 


/\w,r',d' _ jL 
A v,r,d ~ 


4’ 


dE(k w r i d ') 
dk 


-l 


x \M v ' w (k V)r)d , k w y )d ’)[ 


fo[E(k w r / d i )) 

_ 1 - f 0 (E(k v ^ d )) 


(23) 


B 


w,r' ,d' 
v,r,d 


= h 


w,r' ,d' 


dE(k w y,d ') 

dk 


-l 


X I M v ’ w (k v , r ,d,k w y, d ,)\‘ 


f 0 (E(k Vi r >d )) 
Jo(E(k w y, d '))\ 


(24) 


where the non-zero entries of the matrix representing the linear problem are governed by 
the Kronecker symbols f , that are dehned so to enforce energy conservation. 

Eq. 22 has been written for a single scattering mechanism. In order to accommodate 
several scattering mechanisms in our calculations, we do not resort to an approximated 
treatment based on the Matthiessen rule, 4 but instead follow Ref. 23 and notice that Eq. 22 
can be written in the concise matrix notation M^g = G, where M^ is a matrix specific of 
the scattering mechanisms s, g is the unknown vector and G is the vector at the right hand 
side of Eq. 22 and consisting of known quantities. Hence the unknown vector g corresponding 
to several scattering mechanisms can be obtained by solving the linear problem 


- N SC 

^M (s) 

8=1 


a = G 


(25) 


where Eq. 22-Eq. 24 will totally define the entries of the matrix M^ for each scattering 
mechanism. 
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IV. RESULTS AND DISCUSSIONS 


Fig. 3 shows the energy contours of the conduction-band in the first BZ for strained single 
layer M 0 S 2 . In an unstrained material the lowest and the second lowest minimum in the 
conduction band are denoted as K-valley and Q-valley, respectively. The 6 Q-valleys are 
degenerate for unstrained and biaxial strain conditions. With the application of uniaxial 
strain, however, they split into 4 Qa- valleys and 2 Qe-valleys with different effective masses 
and energy minima. Fig. 4 illustrates the bandstructure of unstrained and strained single 
layer M 0 S 2 including K, Qa, and Qb valleys. Under compressive strain one of the Qa or Qb 
valleys becomes the lowest valley. 

The energy distance between these K-valley and Q-valley for unstrained material is eval¬ 
uated to be 160 meV, in agreement with Ref. 10. Tensile strain increases this energy 
distance, which is instead reduced by a compressive strain. In particularly, a relatively 
large compressive strain lowers the energy of Q-valley so that it becomes the lowest valley 
as shown in Fig. 5(a)-(c). Here we can anticipate that, while under tensile strain one can 
neglect the scattering between Q and K-valleys, under compressive strain this type of scat¬ 
tering can significantly affect the mobility. Assuming a non-parabolic dispersion relation 
E( 1 + aE ) = h 2 kf/2m* l + h 2 k^/2m^, the longitudinal m* and transverse m* t effective mass 
and also the non-parabolicity factor a are extracted from the DFT-calculated electronic 
bandstructure and reported in Fig. 5(d)-(f). As can be seen in Fig. 5(a)-(c), under com¬ 
pressive uniaxial strain the energy minima of all K- and Q-valleys are quite close, while at 
large compressive biaxial strain the K-valley lie at higher energy and their contribution to 
mobility can be neglected. 

We compare in Table IV our calculated mobilities at various carrier concentrations with 
the experimental data reported in Ref. 25 for unstrained single-layer M 0 S 2 embedded be¬ 
tween SiCU and HfCU with impirity density 4 xlO 12 . At T = 100 K the effect of piezoelectric 
can be ignored. 11 Very good agreement with experimental data validates the bandstructure 
and mobility models employed in this work. 

The strain-dependency of intrinsic phonon limited mobility is presented in Fig. 6(a). 
Apparently, the effects of compressive and tensile strain on mobility are very different, which 
can be mainly explained by considering the role of inter-valley scattering. For example, 
with tensile strain the minimum energies of Qa and Q^-valleys are much higher than that 
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of K-valley, which suppresses inter-valley scattering. Under compressive strain, instead, the 
inter-valley scattering cannot be neglected because of the smaller energy difference between 
these valleys. With tensile biaxial strain, the mobility increases because of the reduction of 
the effective mass and also the increase of the energy difference between K and Q-vallcys, 
which results in the reduction of the inter-valley scattering rate. With a tensile biaxial 
strain of 5% the phonon limited mobility becomes 75% higher than that of unstrained 
material. In contrast, a compressive biaxial strain of 0.8% strongly reduces the mobility 
due to the reduction of energy difference between K and Q-valleys (see Fig. 5(a)) and 
increased inter-valley scattering. With further increase of compressive biaxial strain, Q- 
valleys become the lowest ones and thus dominate the mobility. At a strain value of about 
2.5% the contribution of K-valleys to mobility becomes negligible and the mobility behavior 
is completely determined by the Q-valleys. Longitudinal and transverse effective masses 
of Q-valleys are not equal and are somewhat changed by strain, however, the different 
angular dependency of mobility along the armchair and zigzag direction tends to compensate 
the changes of effective masses and the overall mobility remains nearly constant at larger 
compressive strain values. 

Under tensile uniaxial strain the mobility is hardly affected by a strain along the zigzag 
direction, while it increases for strain along the armchair direction. In both cases the vari¬ 
ation of the effective mass and non-parabolicity factor with strain determine the mobility 
behavior. Under a compressive uniaxial strain along the armchair direction, Qa becomes 
the lowest valley, while for a strain along the zigzag direction Qb is the lowest one. These 
results emphasize that the contribution of both Qa and Qb valley should be included for an 
accurate calculation of mobility. Under a compressive strain of about 1.5% the mobilities are 
strongly reduced, but they remain nearly constant for larger strain magnitudes. Moreover, 
we notice that for a strain along the zigzag direction, the mobility along the strain direction 
becomes slightly larger than the mobility in the armchair direction. 

Fig. 6(b) reports the mobility in the presence of intrinsic phonon and charged impurity 
scattering. The top and bottom oxide are assumed to be Si0 2 and both carrier and impurity 
concentrations are 10 12 cm -2 . Except for a global reduction of the mobility, the behavior 
of the mobility with strain is similar to Fig. 6(a) corresponding to phonon limited mobility. 
The results presented in Fig. 6(c) correspond to the same parameters as in Fig. 6(b), except 
for a reduction of carrier concentration to 10 11 cm -2 . As the carrier concentration decreases 
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the effect of static screening becomes weaker and the mobility is further reduced. Fig. 6 (d) 
illustrates the mobility as a function of strain with the same parameters used in Fig. 6 (b), 
expect for the top and bottom gate oxide which is AI 2 O 3 . A high-ft dielectric implies a larger 
dielectric screening and increases the mobility. Under this condition, with a tensile biaxial 
strain of 5% and a tensile uniaxial strain of 5% along the armchair direction the mobility 
increases by 53% and 43%, respectively, compared to an unstrained single-layer M 0 S 2 . For a 
better comparison, Fig. 7 shows the room temperature mobility versus carrier concentration 
and also versus the dielectric constant for the unstrained material and for 5% tensile strain 
in either a biaxial or a uniaxial configuration along the armchair direction with an impurity 
density equal to 10 12 cm -2 . As can be seen in Fig. 7(a), because of screening the mobility 
increases with the carrier concentration for both unstrained and strained cases. Fig. 7(b) 
indicates that the strain induced mobility enhancement with high-ft dielectric materials is 
slightly larger than that with low-ft materials. 


The effect of unscreened and screened remote phonon scattering on the mobility of un¬ 
strained and 5% biaxial strained single layer M 0 S 2 are compared in Fig. 8 . Except for a 
global increase of mobility values. The mobility dependence on the dielectric constant ft is 
not significantly affected by the screening of SO phonons. As can be seen, for relatively small 
ft values, mobility improves with increasing ft because of the dielectric screening of charged 
impurities. 1 ' At high ft values, however, the mobility decrease with increasing ft because the 
corresponding smaller SO phonon energies (see Table III) tend to increase momentum re¬ 
laxation time via SO phonons. For the conditions considered in Fig. 8 (temperature, carrier 
and impurity concentrations, and semi-infinite dielectrics with Si 02 as the bottom oxide), 
AIN appears to be the optimal top dielectric material for strained and also unstrained single 
layer M 0 S 2 . Fig. 9 shows the temperature dependency of the mobility for unstrained and 
5% biaxial strain with HfCU as the top oxide. As expected the effect of inelastic remote 
phonons increases with temperature for both unstrained and strained cases. Therefore, it is 
expected that the optimal material as a top dielectric for temperatures above(bcllow) 300 
K, should have a lower(higher)-ft compared to AIN. 
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V. CONCLUSION 


A comprehensive theoretical study on the role of strain on the mobility of single-layer 
MoS 2 is presented. DFT calculations are used to obtain the effective masses and energy 
minima of the contributing valleys. Thereafter, the linearized BTE is solved for evaluating 
the mobility, including the effect of intrinsic phonons, remote phonons, and screened charged 
impurities. The results indicate that, a tensile strain increases the mobility, while compres¬ 
sive strain reduces the mobility. Furthermore, biaxial strain and uniaxial strain along the 
armchair direction increase the mobility more effectively. The strain-dependency of the mo¬ 
bility of MoS 2 is rather complicated and strongly depends on the relative positions of Q and 
K-valleys and the corresponding inter-valley scattering. The presented results pave the way 
for a possible strain engineering of the electronic transport in MoS 2 based electron devices. 
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TABLE I. Deformation potentials for inelastic phonon assisted transitions in single layer M 0 S 2 . 
All parameters are taken from Ref. 16. 


Phonon momentum Electron transition Deformation potential 


r 

K— 

D ac = 4.5 eV 

r 

K— 

D 0p = 5.8 x 10 8 eV/cm 

K 

K—^K 7 

D ac = 1.4 x 10 s eV/cm 

K 

K—^K 7 

D op = 2.0 x 10 s eV/cm 

Q 

K^Q 

D ac = 9.3 x 10 7 eV/cm 

Q 

K^Q 

D op = 1.9 x 10 8 eV/cm 

M 

K—>• Q 

D ac = 4.4 x 10 7 eV/cm 

M 

K-» Q 

D op = 5.6 x 10 8 eV/cm 


r 

r 

Q 

Q 

M 

M 

K 

K 

Q 

Q 

M 

M 


Q -> Q 
Q -> Q 
Q ->Q 
Q ->Q 
Q 

Q ->Q 
Q -> Q 
Q -> Q 
Q ->K or K' 
Q ->K or K' 
Q —)-K or K' 
Q —)-K or K' 


D ac = 2.8 eV 


D op = 

= 7.1 

X 

10 8 

eV/cm 

D ac = 

= 2.1 

X 

10 s 

eV/cm 

D op = 

= 4.8 

X 

10 8 

eV/cm 

Dac = 

= 2.0 

X 

10 s 

eV/cm 

Dop = 

= 4.0 

X 

10 8 

eV/cm 

Dac = 

= 4.8 

X 

10 8 

eV/cm 

Dop = 

= 6.5 

X 

10 8 

eV/cm 

Dac = 

= 1.5 

X 

10 s 

eV/cm 

Dop = 

= 2.4 

X 

10 8 

eV/cm 

D ac = 

= 4.4 

X 

10 s 

eV/cm 

Dop = 

= 6.6 

X 

10 8 

eV/cm 
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TABLE II. Phonon energy for intra-valley and inter-valley transitions at the K, M, and Q points 
of single layer M 0 S 2 as reported in Ref. 16. As discussed in Ref. 16, the energy values for acoustic 
(optical) phonon modes is the average of phonon energies of transverse and longitudinal (transverse, 
longitudinal and homo-polar) modes. 


Phonon mode 

r 

K 

M 

Q 

Acoustic phonon energy [meV] 

0 

26.1 

24.2 

20.7 

Optical phonon energy [meV] 

49.5 

46.8 

47.5 

48.1 
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TABLE III. Parameters for the dielectric materials taken from (a) Ref. 12 and (b) Ref. 24 and 
corresponding calculated SO phonon frequencies hco SOt t ox and huj so ^ ox . In all of the cases an Si02 
bottom oxide is assumed. 


Top oxide dielectric material 

Si0^ a) 

BN^ 

A1N(°) 

Al 2 oj“ } 

Hf0^ a) 

Zr0^ a) 

fc tox 

3.9 

5.09 

9.14 

12.53 

23 

24 

oo 

fc tox 

2.5 

4.1 

4.8 

3.2 

5.03 

4 

WTO, tox [meV] 

55.6 

93.07 

81.4 

48.18 

12.4 

16.67 

w so ,tox [meV] (Evaluated in this work) 

69.4 

100.5 

104.3 

83.9 

21.3 

30.5 

w so ,box [meV] (Evaluated in this work) 

69.4 

60.1 

58.0 

54.2 

61.1 

62.9 
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TABLE IV. Comparison of the calculated mobility in this work with the experimental data of 
Ref. 25. T = 100 K and the impurity density is 4 x 10 12 cm -2 . 


Carrier concentration [cm -2 ] 

7.6 x 10 12 

9.6 x 10 12 

1.15 x 10 13 

1.35 x 10 13 

Calculated mobility, this work [cm 2 /(Vs)] 

93 

106 

114 

122 

Experimental mobility [cm 2 /(Vs)] 

96±3 

111±3 

128±3 

132±3 
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LIST OF FIGURES 


Fig. 1: (a) Equi-energy contours in the first Brillouin zone for the unstrained single layer 
MoS 2 . The angle 9 that describes the Qa valleys orientation in /c-space is also depicted 
in the figure. It should be recalled that the zigzag direction in fc-space corresponds to 
the armchair direction in real space, (b) The bandstructure of unstrained single layer 
MoS 2 in the first Brillouin zone and along the symmetry directions that are illustrated 
in (a). 

Fig. 2 : Illustration of several phonon assisted inter-valley transitions in single layer MoS 2 
for (a) transitions from K-valley to other valleys; (b) transitions from Q^-valley to 
other valleys; (c) transitions from Q^-valley to other valleys. The figure also sets the 
notation used in Table I and Table II to identify phonon assisted transitions. 

Fig. 3: Equi-energy contours for single layer MoS 2 under: (a) compressive biaxial strain; 

(b) tensile biaxial strain; (c) compressive uniaxial strain along the armchair direction; 
(d) tensile uniaxial strain along the armchair direction; (e) compressive uniaxial strain 
along the zigzag direction; (f) tensile uniaxial strain along the zigzag direction, (g) 
Extracted effective mass of K-valley along all directions in polar coordinate for un¬ 
strained MoS 2 and under tensile biaxial and uniaxial strain along armchair and zigzag 
directions. The nearly circular shape of the effective mass plot justifies the assumption 
of isotropic bandstructure. The strain magnitude is 4% in all strained cases. The lon¬ 
gitudinal and transverse effective masses of Q-valleys vary with the strain conditions. 

Fig. 4 : The band structure of unstrained and strained single layer MoS 2 . BI: biaxial strain, 
UA: uniaxial strain along armchair direction; UZ: uniaxial strain along zigzag direction. 
The strain magnitude is 4% in all strained cases. 

Fig. 5 : The minimum energies of valleys (solid-lines) and the angle 9 (dotted lines) between 
the longitudinal direction of Qa valleys and zigzag direction in Uspace as illustrated 
in Fig. 1(a) under: (a) biaxial strain; (b) uniaxial strain along the armchair direction; 

(c) uniaxial strain along the zigzag direction. The 9 angle in Fig. 5(a)-(c) corresponds 
to the Qa valley indicated in Fig. 1(a), and the 9 angle of the other Qa valleys can 
be inferred from symmetry considerations. The 9 angle for Q^ valleys has a negligible 
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dependence on strain (not shown) and it is approximately zero (see Fig. 1(a)). The 
effective masses (solid-lines for longitudinal and dashed-lines for transverse) and the 
non-parabolicity factor (a) (dotted-lines) of various valleys under: (d) biaxial strain; 
(e) uniaxial strain along the armchair direction; (f) uniaxial strain along the zigzag 
direction. The longitudinal and transverse effective masses of K-valley are assumed to 
be equal. 

Fig. 6: (a) Phonon limited mobility of single layer M 0 S 2 as a function of strain with a 
carrier concentration n = 10 12 cm -2 . Mobility limited by phonon and screened charged 
impurity scattering with SiC >2 as the gate oxide (e r = 3.9) and carrier (n) and charged 
impurity concentration (ui mp ) for: (b) n = n; mp = 10 12 cm -2 ; (c) n = 10 11 cm -2 and 
77-imp = 10 12 cm -2 , (d) Same as (b), except for the gate oxide which is A1 2 0 3 . In 
the legend, BI, UA, and UZ denote biaxial strain, uniaxial strain along the armchair 
direction, and uniaxial strain along the zigzag direction respectively. The subscripts A 
and Z indicate the component of the mobility along the armchair or zigzag direction. 
For example: UZ^ is the mobility along armchair direction for a uniaxial strain along 
zigzag direction. 

Fig. 7: (a) The mobility versus carrier concentration with and without screening for the 
unstrained M 0 S 2 , for a tensile biaxial strain of 5%, and for a uniaxial strain of 5% 
along the armchair direction. ni mp = 10 12 cnr 2 . (b) The mobility versus the relative 
dielectric constant for unstrained M 0 S 2 and for strain conditions as in (a), n — ni mp = 
10 12 cnr 2 . The strain induced mobility enhancement is shown on the right-side of the 
y- axis. 

Fig. 8 : The mobility accounting for intrinsic phonon and charged impurity scattering (tri¬ 
angle), and for either unscreened (rectangle) or screened (circle) SO phonon scattering 
as a function of top oxide dielectric constant for unstrained (blue line) and 5% biax¬ 
ial strain (red line). Numbers 1 to 6 indicate the k value corresponding to dielectric 
materials studied in this work (see also Table III). I 11 particular, (1): Si 02 , (2): BN, 
(3): AIN, (4): AI 2 O 3 , (5): Hf 02 , and ( 6 ): Zr 02 - In all cases the back oxide is assumed 
to be Si 02 - T = 300 K, the impurity and carrier concentrations are equal to 4 x 10 12 
cnr 2 and 10 13 cnr 2 , respectively. These values are consistent with experimental data 
reported in Ref. 25. 
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Fig. 9: The mobility with the inclusion of intrinsic phonon and charged impurity scattering 
(dash line) and with the inclusion of screened SO phonon (solid line) versus tempera¬ 
ture for a Si0 2 /MoS 2 /Hf0 2 structure for unstrained (blue line) and 5% biaxial strain 
(red line). The impurity and carrier concentrations are equal to 4 x 10 12 cm -2 and 
10 13 cm -2 , respectively. 
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